%FIR滤波器
%omega_p=0.6pi omega_s=0.4pi, ap<=0.3db , as>=40db 高通
%使用一型，因为四种都可以实现，一般根据阻带衰减选择窗，没有给出，任选一个。
clear; clc;
% win=rectwin(L);
% win=hamming(L);
wp=0.6*pi;
ws=0.4*pi;
ap=0.3;
as=40;

wc=(wp+ws)/2;
N=ceil(6.2/(wp-ws)); %算滤波器长度,根据汉宁窗过渡带宽度反求N
N=mod(N+1,2)+N;
m=N-1;
k=0:m;
w=hanning(N)';
hd=-(wc/pi)*sinc(wc*(k-0.5*m)/pi);
hd(0.5*m+1)=hd(0.5*m+1)+1;%转换为高通
h=hd.*w;
omega=linspace(0,pi,512);
mag=freqz(h,[1],omega);%频率响应
plot(omega/pi,20*log10(abs(mag)))


